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Abstract 

Modifications to polar-gap models for pulsars are discussed for the case where the 
surface magnetic field, B s , of the neutron star is strong. For B > 4 x 10 8 T, the curvature 
7-quanta emitted tangentially to the curved force lines of the magnetic field are captured 
near the threshold of bound pair creation and are channelled along the magnetic field as 
bound electron-positron pairs (positronium). The stability of such bound pairs against 
ionization by the parallel electric field, in the polar cap, and against photoionization 
is discussed. Unlike free pairs, bound pairs do not screen En near the neutron star. As 
a consequence, the energy flux in highly relativistic particles and high-frequency (X-ray 
Q-f and/or 7-ray) radiation from the polar gaps can be much greater than in the absence of 
positronium formation. We discuss this enhancement for (a) Arons-type models, in which 
particles flow freely from the surface, and find any enhancement to be modest, and (b) 
Ruderman-Sutherland-type models, in which particles are tightly bound to the surface, 
and find that the enhancement can be substantial. In the latter case we argue for a self- 
consistent, time-independent model in which partial screening of E\\ maintains it close to 



$_i ■ the threshold value for field ionization of the bound pairs, and in which a reverse flux 



of accelerated particles maintains the polar cap at a temperature such that thermionic 
emission supplies the particles needed for this screening. This model applies only in a 
restricted range of periods, P2 < P < Pi, and it implies an energy flux in high-energy 
particles that can correspond to a substantial fraction of the spin-down power of the 
pulsar. 

Nonthermal, high-frequency radiation has been observed from six radio pulsars and 
Geminga is usually included as a seventh case. The nonthermal luminosity can be higher 
than can be explained in terms of conventional polar-gap and outer-gap models. The self- 
consistent polar-gap model proposed here alleviates this difficulty, provided the magnetic 
field satisfies B >4x 10 8 T (which is so for five of these pulsars, and plausibly for the other 
two if a modest nondipolar component is assumed), and the surface temperature (in the 
absence of heating by the reverse flux) satisfies T s < 10 6 K, so that thermionic emission 
from the surface is unimportant. It is argued that sufficient power is available to explain 
the observed high-frequency radiation of most of these pulsars. However, the Crab and 
PSR 0540-69 have periods P < P2, and we suggest that an outer-gap model is more 
appropriate for these. 
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1. INTRODUCTION 



The creation of electron-positron pairs by decay of 7 rays as they propagate across magnetic 
field lines is an essential ingredient in the population of a pulsar magnetosphere with 
plasma. In the absence of plasma, a strong electric fields E results from the rotation 
of the magnetized neutron star (Goldreich and Julian 1969). The vacuum field has a 
nonzero parallel component, E\\ = (E ■ B)/|B| along the magnetic field B, and this E\\ can 
accelerate primary particles to ultrarelativistic energies (e.g., the review by Michel 1991). 
The source of the primary particles is different in different models. In the absence of free 
ejection from the stellar surface, the primary particles come from a pair cascade initiated, 
for example, by decay of stray 7 rays into pairs (Ruderman and Sutherland 1975). In 
such models the vacuum E\\ is present immediately above the surface and is said to form 
a vacuum gap. In models with free ejection of particles from the stellar surface, these 
particles screen E\\ immediately above the surface, and set up the corotation electric field 
E ro t = — (O x r) x B, where O is the angular velocity of rotation. The divergence of E rot 
requires a charge density, en GJ , where n GJ is the Goldreich- Julian density (Goldreich and 
Julian 1969). In such models an E\\ develops, forming a vacuum-gap-like region above the 
surface, due to the actual charge density deviating increasingly with increasing height from 
the Goldreich- Julian value, n 7^ n GJ . In either model, the primary particles emit 7 rays, 
due to curvature emission and other processes. Some of these 7 rays are absorbed in the 
magnetic field by creating free electron-positron pairs (7 + B — > e + + e~ + B). Provided 
these secondary particles are created in the gap, the E\\ field can accelerate the electron 
and the positron in opposite directions, allowing a net charge density to build up, and 
it is this charge density of the secondary particles that screens E\\ . Regions where E\\ is 
unscreened are called gaps, and a gap that forms near the magnetic poles of the pulsar 
is called a polar gap. The height, H, of a polar gap is defined by the height above the 
stellar surface where the charge density due to the charge-separated pairs is sufficient to 
screen E\\. Besides polar-gap models there are also outer-gap models (e.g., Cheng, Ho and 
Ruderman 1986a,b), where the region with £j| 7^ occurs far from the stellar surface, 
where the magnetic field is much weaker than in the polar-cap regions. In comparing the 
slot-gap with polar-gap and outer-gap models, we note that the slot gaps form on the 
boundary of the closed S-field region (Arons 1983). Primary particles which are ejected 
from the neutron star surface and are accelerated in the slot gap gain the main part of 
their energy far from the neutron star. In this respect, slot gaps differ significantly from 
the polar gaps and may be regarded as intermediate between polar gaps and outer gaps. 

Polar gaps are thought to be a source of the energy for nonthermal radiation of pulsars 
(Sturrock 1971; Ruderman and Sutherland 1975; Michel 1975; Arons 1979; Arons and 
Scharlemann 1979; Cheng and Ruderman 1980; Arons 1981; Mestel 1981, 1993; Mestel et 
al. 1985; Fitzpatrick and Mestel 1988a,b; Shibata 1991). In particular, in most models the 
radio emission is attributed to processes in polar gaps. However, the radio luminosities of 
the pulsars are small compared with the total power loss estimated from the observed spin 
down and associated rotational energy loss ( < 10" 5 ). Most of the power loss is attributed 
to relativistic particles created as secondary pairs in the polar gaps and escaping in a 
relativistic wind (Ruderman and Sutherland 1975; Arons 1979; Michel 1991). In most 
pulsars the only nonthermal radiation observed is in the radio range and, apart from the 
slowing down, the large inferred power loss has no direct observational signature. However, 



there are a few radio pulsars that do have high-frequency emission (in X rays and/or 7 
rays, we do not distinguish between them in this paper) with a nonthermal luminosity that 
is a substantial fraction of the inferred rotational power loss. Nonthermal high-frequency 
radiation requires high energy particles, and it is widely assumed that the high-frequency 
luminosity correlates with the power in primary particles (e.g., Harding and Daugherty 
1993). The power in the primary particles, Ar pr i m eA<^, is limited by the rate of injection 
of such particles, A^ pr i m , which is less than n GJ c times the area of the polar cap, and 
the potential energy, eA<p ~ eE^H, that particles gain in crossing the polar gap. The 
maximum potential difference across the polar gap is when there is no screening at all: 

AVW = - (6.6 x 10 14 V) ( -§-) (t^Y\ (1-1) 
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where is the dipole component of the magnetic field at the magnetic pole on the neutron 
star surface, O = 2tv/P is the angular speed of rotation of the neutron star and R ~ 10 4 m 
is the radius of the neutron star. For Aip = A</? max the implied maximum value for power 
in primary particle is equal (within a factor of order unity) to the rotational energy loss. 
For existing polar-gap models the value of A(p is, as a rule, considerably less than A</? max . 
In Arons-type models, in which particles flow freely from the surface, this is due both (a) to 
E\\ being much less than the vacuum value due to screening by the primary particles that 
flow freely from the surface, and (b) to the height being less than the maximum value, 
H ~ R, for such a model due to screening by secondary pairs. In Ruderman- Sutherland- 
type models, in which particles are tightly bound to the surface, A(p <C A</? max results 
from H <C AR P due to screening by secondary pairs, where AR P is the radius of the polar 
cap which is the maximum value of H in this kind of polar-gap model. 

Strong high-frequency radiation is observed from the galactic pulsars PSR 0531+21, 
PSR 0833-45, PSR 1055-52, PSR 1509-58 and PSR 1706-44, and PSR 0540-69 in the 
Large Magellanic Cloud (e.g., the review by Ulmer 1994). The observed radiated power 
for all these pulsars is concentrated in the X-ray or 7-ray ranges. The 7-ray pulsar Geminga 
is probably also a radio pulsar (Halpern and Holt 1992) which is "radio quiet" because 
its radio beam does not intersect the Earth (Ozernoy and Usov 1977). The high 7-ray 
luminosities are inferred assuming emission into a large solid angle, ~ 2-7rsterad (Ulmer 
1994). For some of the the pulsars, notably for PSR 1055-52, among the present polar-gap 
models the only one that appears capable of explaining the inferred very high luminosities 
is that of Sturrock (1971). In the Sturrock model the potential drop has the maximum 
possible value, and the nonthermal power is of order the rotational power loss. However, 
the Sturrock model is not self-consistent because it neglects the finite limit on the height 
of the polar gap due to screening. In the Sturrock model the free electron-positron pairs 
created by the curvature 7-quanta inside the polar gap are assumed not to feel the electric 
field E\\ and not to screen this field so that the full potential drop A</? max is implicitly 
available. This internal inconsistency implies that the Sturrock model should not be used 
for detailed estimates of the pulsar luminosities. 

Ruderman and Sutherland (1975) were the first to develop a self-consistent polar-gap 
model in which the screening of the electric field E\\ by the electron-positron pairs created 
in it is taken into account. Consideration of this screening led Ruderman and Sutherland 
(1975) to conclude that the potential across the polar gap cannot exceed A(p RS ~ a few 
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xlO V. This upper limit on the potential across the polar gap is valid for any polar-gap 
model in which free pairs are created by 7-quanta absorption in the magnetic field. For 
young, rapidly-rotating pulsars with typical parameters, P < 0.1 s and Bg > 10 8 T, one has 
A<p RS < 10" 3 A^ max . The total power in primary particles is proportional to Ay?, and the 
maximum luminosity expected in all existing polar-gap models (Ruderman and Sutherland 
1975; Arons 1979, 1981; Arons and Scharlemann 1979; Cheng and Ruderman 1980; Mestel 
et al. 1985) is smaller than the inferred nonthermal luminosity in X-rays and 7 rays (section 
5). 

In all existing polar-gap models, the formation of pairs in the pulsar magnetospheres is 
due to the single-photon mechanism, 7 + 5 — > e + +e~ +B, and the resulting pairs are free. 
However, this assumption that the pairs are free is not valid if the magnetic field is strong 
enough, specifically for B > 0.1 5 cr , where B cr = m 2 c 2 /eh = 4.4 x 10 9 T is known as the 
critical field. In such a strong magnetic field, the curvature 7-quanta emitted tangentially 
to the curved force lines of the magnetic field are captured near the threshold of bound 
pair creation and are then channelled along the magnetic field as positronium, that is, as 
bound pairs (Shabad and Usov 1982, 1985, 1986; Herold, Ruder and Wunner 1985; Usov 
and Shabad 1985; Meszaros 1992; Shabad 1992 and references therein). This positronium 
may be stable in the polar gaps against both the ionizing action of the electric field and 
against photo-ionization (Shabad and Usov 1985; Bhatia, Chopra and Panchapakesan 
1988, 1992). Unlike free pairs, such bound pairs do not screen the electric field E\\ near 
the pulsar. Screening requires a net charge density, which can build up due to free pairs 
being separated by £j|, but cannot build up if the pairs remain bound. As a result the 
height of the gap, which is determined by the height at which screening becomes important, 
is greater than it would be in the absence of formation of positronium. The assumption 
that the power in primary particles is proportional to H in polar-gap models, implies that 
this total luminosity increases when positronium formation becomes important (Usov and 
Shabad 1985). We argue in section 5 that the nonthermal luminosity may even become 
comparable with the maximum possible, as in the model of Sturrock (1971). 

Our primary objective in this paper is to discuss how models of pulsar electrodynamics 
need to be modified when the magnetic field is strong enough for creation of bound pairs 
to dominate over the formation of free pairs. First, however, we review several important 
preliminary aspects of the problem. The rate of injection of the primary particles depends 
on the properties of the matter at the surface of a magnetic neutron star. The particle 
outflow from the neutron star surface is discussed in section 2. Particle acceleration and 
generation of curvature 7 rays in the pulsar magnetospheres are described in section 3. 
Propagation of 7 rays, including photon splitting, and creation of electron-positron pairs 
in a strong magnetic field are discussed in section 4. In section 5, the polar-gap model 
is developed for pulsars with strong magnetic fields at their surface, B s > 0.15 cr , and 
the maximum value of the nonthermal luminosity is estimated for such a pulsar. The 
interpretation of the observational data and some theoretical predictions on the nonthermal 
radiation of pulsars with B s > 0.1 B CT are also given in section 5. A discussion and 
summary are presented in section 6. 

2. THE NEUTRON STAR SURFACE AND PARTICLE EJECTION 

The electric field distribution and particle acceleration in the polar gaps of pulsars depend 
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on the character of particle outflow from the pulsar surface. In some models the binding to 
the surface is assumed sufficiently strong that no particles are ejected from the surface, and 
in other models the binding is assumed sufficiently weak that particles (either electrons 
or ions) flow freely from the stellar surface. Particle ejection from the surface depends 
on the composition of the neutron star surface, the structure of the surface with a strong 
magnetic field and the surface temperature T s . In this section, after a brief review of 
polar-gap models we discuss the properties stellar surface and their effect on the binding 
or ejection of ions. We then discuss models involving free ejection of electrons and free 
ejection of ions. 

2.1. Polar-gap models 

There are several kinds of polar-gap model. These may be classified in two ways: whether 
ions or electrons tend to be ejected from the surface, and whether E\\ is zero or nonzero 
at the stellar surface. 

The sign of the charge of the particles that tend to be ejected from the neutron star 
surface depends on the sign of O • B. Electrons tend to be ejected for O • B > and 
ions for ft ■ B < 0. The argument for this (Goldreich and Julian 1969; Ruder man and 
Sutherland 1975; Arons 1979, 1981) is that plasma tends to corotate with the star, and 
the divergence of the corotation electric field implies a charge density, en GJ , where 

= ^-< 7xl ° 17m " 3 »(^)(o^)" lcose ' (21) 

is the Goldreich-Julian density, and is the angle between the angular velocity ft and the 
magnetic field B. Thus the tendency to corotate requires positive charges above the polar 
caps for ft • B < and negative charges over the polar caps for ft ■ B > 0. The required 
charge density may be produced in two ways: ejection of particles with the appropriate 
sign of the charge from the stellar surface, or acceleration (by E\\) from the magnetosphere 
back to the stellar surface of charges of the opposite sign. 

The most familiar model in which there is no ejection of particles from the stellar 
surface is that of Ruderman and Sutherland (1975). This model applies only to neutron 
stars with ft ■ B < (neutron stars with ft ■ B > are assumed not to be pulsars). The 
absence of ion flow from the surface is attributed to an assumed strong binding of iron 
nuclei. In the Ruderman-Sutherland model, the field E\\ is maximum at the surface and 
decreases with distance. All the particles in the pulsar magnetosphere are created from pair 
production through 7-ray absorption. The magnetosphere is populated by a pair plasma 
created near the polar cap, and the (positive) Goldreich-Julian charge density is created by 
electrons being accelerated back to the stellar surface, leaving a net excess (and outward 
flux) of positrons. A model with somewhat similar properties is that of Beskin, Gurevich 
and Istomin (1986) in which the pair plasma is assumed to be created in a postulated 
double layer immediately above the surface on the star. 

In contrast, in the model of Arons (1979, 1981) it is assumed that charged particles 
(electrons or ions depending on the sign of ft • B) flow freely from the neutron star surface. 
In this model the electric field E\\ is equal to zero at the surface, and increases with 
distance above the surface. In the absence of any additional source of charge, conservation 
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of the current implies a current density J oc B. For so-called favorably curved field lines, 
the ratio of the charge density (n = J/ec oc B) to the Goldreich- Julian charge density 
(n GJ oc ft • B/fi) decreases with distance above the stellar surface. This implies that the 
screening of E\\ is incomplete, and E\\ increases toward the vacuum value. The resulting E\\, 
which is smaller than that in the Ruderman- Sutherland model, causes the ejected electrons 
or ions to become highly relativistic and causes a net flux of pair-produced electrons or 
positrons back to the stellar surface to allow the charge density to be maintained at the 
Goldreich-Julian value. Mestel et al. (1985) pointed out that for the unfavorably curved 
field lines, along which n/n GJ increases, an E\\ of opposite sign to the vacuum field develops 
and plays essentially the same role as the E\\ in the Arons model. A third kind of polar-gap 
model (Cheng and Ruderman 1980) is an intermediate case where the particles flow from 
the pulsar surface but not freely. In such a model the field E\\ is nonzero at the pulsar 
surface but it is smaller than in the model of Ruderman and Sutherland (1975). 

In the discussion below, we use the Ruderman- Sutherland and Arons models for illus- 
trative purposes. 

2.2. Composition of the neutron star surface 

The strength of the binding of ions to the stellar surface is an important ingredient in 
pulsar models with f2 ■ B < 0, and this binding depends on the ionic composition. We 
note the following arguments for and against the composition being primarily Fe 56 . 

1) . Rosen and Cameron (1972) showed that the outermost layers of the neutron star 
atmosphere, in which no mass is being ejected, are composed almost entirely of He with 
trace amounts of Fe 56 . The He 4 is produced by photo dissociation of the iron-peak nuclei 
(assumed to be the major initial constituent) at the high temperatures encountered during 
collapse to the neutron star state. The total mass of the He 4 at the neutron star surface 
is AMn e ~ 1.3 x 10 18 kg ~ 10 -12 Mq. However, the luminosity of a very young neutron 
star is highly super-Eddington, and so radiation pressure drives matter away from the 
neutron star surface. After the loss of about 10 _12 Mq, only Fe 56 remains at the surface 
(Rosen and Cameron 1972). For young neutron stars the mass-loss rate is as high as 
~ 0.005 M©s _1 , and the mass ejected from the neutron star surface is of order 10~ 3 - 
10 _2 M Q (Woosley and Baron 1992; Levinson and Eichler 1993). Thus the ejected mass is 
3> AMn e , suggesting that the neutron star surface consists almost entirely of Fe 56 . 

2) . The composition of the polar caps of pulsars may differ substantially from the compo- 
sition of the main part of the neutron star surface due to the bombardment by energetic 
particles from the magnetosphere. Pair particles created near the top edge of the polar 
gap and accelerated by the E\\ field back to the stellar surface form showers as they cross 
the surface (Cheng and Ruderman 1977; Jones 1978, 1979). Protons and spallation nuclei 
are produced at the polar caps of pulsars by hadronic photoabsorption of shower photons. 
The flux of reversed electrons (or positrons) through unit surface area of the polar cap is 

F r = £n GJ c, (2.2) 

where £ is a dimensionless parameter, with £ ~ 1 for the model of Ruderman and Suther- 
land (1975) and £ ~ (QR/c) ~ 2 x 10" 3 (P/0.1 s)" 1 < 1 for the model of Arons (1979, 
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1981). The energy of reversed particles is Y r mc 2 with Y r ~ 10 7 . Shower 7 rays pene- 
trate into the neutron star surface layers to ~ 30 radiation lengths (Jones 1978; Bogovalov 
and Kotov 1989) and break up atomic nuclei. In the case of Fe 56 , for which the radi- 
ation length is l r ~ 140kgm~ 2 , 30 radiation lengths corresponds to a column density 
a = 30/ r ~ 4.2 x 10 3 kgm -2 . The mean number of photoabsorption events per one re- 
versed particle with energy ~ 10 7 mc 2 is ~ 10 3 (Hayward 1965; Jones 1978). Hence, from 
equations (2.1) and (2.2), the characteristic time for spallation of all the Fe 56 in the surface 
layer, r = (a / Am p ) / 10 3 F r with A = 56, reduces to 



which is very much shorter than the age of any pulsar. This suggests that all the Fe 56 
should be destroyed in the surface layers of the polar caps. 

3). There are several contrary arguments to the foregoing, (a) The most important 
photoabsorption reaction is the formation of the giant dipole resonance (e.g., Hayward 
1965; Jones 1978), which, for iron-group nuclei, decays predominantly through neutron 
emission. Most of the photons with energies of ~ 15 - 30MeV mainly responsible for 
spallation of Fe 56 are generated in the deep layer, > 10 Z r , (e.g., Jones 1978; Bogovalov 
and Kotov 1989). This suggests that the surface layer itself is relatively unaffected by the 
bombardment, (b) The resulting neutrons are distributed more or less uniformly inside the 
surface layer due to diffusion (Jones 1978), and are recaptured by the other nuclei. Neutron 
capture can reform Fe 56 , so that it is not permanently destroyed, (c) If the magnetic field 
is strong enough, B s > 5 x 10 8 T, so that solid magnetic metal is formed at the pulsar 
surface (see below), both the generation of photons with energies of ~ 15 - 30MeV and 
spallation of Fe 56 may be suppressed over a few l r due to the Landau- Pomeranchuk effect 
(Rozental and Usov 1985), that is, due to the suppression of the formation of softer photons 
and pairs in the strong magnetic field, (d) For pulsars with f2 ■ B < at the polar caps, 
so that ions are the particles ejected from the stellar surface, ejection of ions can dominate 
over spallation. If so, the spallation products are removed as fast as they are produces, 
and the surfaces of the polar caps consist mainly of Fe 56 . 

In summary, the composition of the surface of the neutron star at its polar caps is 
uncertain. Below we let the mass number, A, of the surface ions be a free parameter, and 
for numerical estimates we assume it corresponds to Fe 56 . 

2. 3. The surface structure of neutron stars with strong magnetic fields 

Whether an Arons-type model or a Ruderman-Sutherland-type model is relevant for the 
polar-cap region of the pulsar magnetosphere depends on whether or not charged particles 
can escape freely from the surface of the neutron star due to thermionic emission. This in 
turn depends on (a) the binding energy, referred to here as the cohesive energy for ions and 
as the work function for electrons, and (b) the surface temperature. The cohesive energy, 
Ae c , was overestimated in earlier literature. 

The structure of matter in the surface layers of neutron stars with B s ^> a 2 B cr ~ 
2.35 x 10 5 T, where a = e 2 /4T7£ fic = 1/137 is the fine structure constant, is largely 
determined by the magnetic field (Ruderman 1971; cf. also Fushiki, Gudmundsson and 
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Pethick 1989 and references therein). It has been suggested that the surfaces of magnetic 
neutron stars with B s ~ 10 8 T consist of a phase of matter where atoms form chains 
aligned along the field lines (Cheng, Ruderman and Sutherland 1974; Flowers et al. 1977). 
The density of the magnetic metal phase is 

^(4x 10 6 kg m" 3 ) ( * V 7 74 )(£) " 3/5 . (2.4) 



10 8 TJ V 56 / V 26 / 

The cohesive energy of the condensed Fe 56 matter was estimated by Flowers et al. (1977), 
using a variational method, to be 
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The properties of both bulk matter and isolated atoms in strong magnetic fields have 
been reconsidered using various different methods, including variational methods (Miiller 
1984; Skjervold and Ostgaard 1984b), the Thomas- Fermi and related methods (Skjervold 
and Ostgaard 1984a), density functional methods (Jones 1985, 1986), the Hartree-Fork 
method (Neuhauser, Langanke and Koonin 1986; Neuhauser, Koonin and Langanke 1987), 
and the Thomas-Fermi-Dirac method with the Weizsacker gradient correction (Abrahams 
and Shapiro 1991). These investigations imply that Flowers et al. (1977) overestimated 
the cohesive energy by a factor of at least a few. The most recent Thomas-Fermi-Dirac 
calculations with the Weizsacker coefficient A = 1/9 (the TFD-^W method) gives cohesive 
energy for Fe 56 to be Ae c ~ 0.91 keV for B = 10 8 T, Ae c ~ 2.9 keV for B = 5 x 10 8 T, 
and Ae c ~ 4.9 keV for B = 10 9 T (Abrahams and Shapiro 1991). The results of these 
calculations may be approximated by 

f B \ °- 73 
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The errors introduced by the TFD-|W method are near one per cent of the binding energy 
of isolated atoms, which is (Abrahams and Shapiro 1991) 

/ Z \ 9/5 , B n 2/5 

£ ^-(55keV)(-) (^j . (2.7) 

From (2.6) and (2.7), it follows that at B ~ 10 s T the errors of the calculations, ~ 10 _2 £;,, 
are of the same order as the cohesive energy. Thus, the qualitative result of Flowers et 
al. (1977) that chains are energetically favored over individual atoms is questionable. It is 
worth noting that the value of A = 1/9 is the standard gradient expansion result for the 
nonmagnetic case (Lieb 1981 and references therein). If the value of A for strong magnetic 
fields is smaller than 1/9, the cohesive energy is smaller than implied by (2.6). Moreover, 
if A is small enough, free atoms are preferred over chains for Z = 26 at B < a few x 10 8 T 
(Neuhauser et al. 1986; Neuhauser et al. 1987). Hence, at B s < a few xlO 8 T there are two 
possibilities: either Fe 56 does not form a magnetic metal at the neutron star surface, or the 
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surface is solid with a small cohesive energy per Fe atom, Ae c < 2 - 3keV (Neuhauser et 
al. 1986). The existence of a magnetic metal is then unimportant for ejection of particles 
from the surface, as discussed below. 

For B > (0.5 - 1) x 10 9 T the surface of a neutron star is probably a magnetic metal, 
provided that the surface temperature is smaller than the melting temperature, T m . An 
estimate of T m is (Slattery, Doolen and DeWitt 1980; Shapiro and Teukolsky 1983; Ogata 
and Ichimaru 1990) 

rm ^ (5 x XO^K) (f ) 2 (A)-" 3 (_|_) 1/3 , (, 8) 

where p, the density of the magnetic metal, may differ from p s , as given by (2.4), by a 
factor two or so (Abrahams and Shapiro 1991). The cohesive energy in such a strong 
magnetic field may be more than (3 - 5) keV. 

2.4- Electron ejection (O ■ B > 0) 

The flow of electrons away from the surface of a pulsar with Q ■ B > 0, assuming the 
neutron star surface to be a magnetic metal, is determined either by thermionic emission 
or by field emission. 

The important condition is that the current density, Jth, from thermionic emission 
provide the Goldreich-Julian charge density. The particles are quickly accelerated to rela- 
tivists energies above the surface, and hence this condition becomes 

Jth > n GJ ec, (2.9) 

Screening due to this charge density then implies E\\ = on the stellar surface. 

For electrons, the current density due to thermionic emission is determined by the 
Richardson-Dushman equation (e.g., Gopal 1974): 

'"^wM-ssO- (2 - 10) 

where m is the electron mass, k is the Boltzmann constant and w is the work function of 
electrons. Typically, the value of w is near the Fermi energy, which is given for a magnetic 
metal by (Ruderman 1971; Flowers et al. 1977) 
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From (2.4) and (2.11), we have 
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Taking w ~ e F and using (2.2), (2.10) and (2.12), we can rewrite the condition (2.9) in 
terms of the characteristic temperature, T e , identified as 

w ( 7 \ 4/5 / B \ 2/5 

r e = 0.04 I .(3.7xlO'K)(-) (^) . (2.13) 

Then setting T = T s as the surface temperature, (2.9) requires T s > T e for thermionic 
emission of electrons to be adequate to provide the Goldreich-Julian density. 

Field emission is relevant only if thermionic emission is inadequate, that is, for T s <T e . 
Then E\\ is nonzero on the surface at the polar caps and the ejection of electrons results 
from quantum mechanical tunneling through the barrier provided by the work function 
and £j|. The current density due to this tunneling is (Beskin 1982) J = ME\\ exp(—N/E\\), 
M = 3 x 10 16 (S/10 8 T)(w/lkeV)- 1 / 2 s- 1 , N = 2 x 10 14 (w/l keV) 3 / 2 V m" 1 . Equating 
this to the current implied by an electron density n GJ outflowing at the speed of light 
allows one to define a characteristic electric field: 

iJe^GxlO^Vm-^— ) , (2.14) 

such that one requires E\\ > E e for field emission to be effective. For E\\ E e , the density 
of electrons ejected due to field emission is negligible. For E\\ > E e the electrons escape 
freely until their density of electrons builds up to ~ n GJ , and they then screen the electric 
field, restricting it to E\\ ~ E e . 

2.5 Ion ejection (O • B < 0) 

The flow of ions away from the solid surface of a pulsar with O • B < is limited by the rate 
of thermionic emission of ions. The density of outflowing ions is (Cheng and Ruderman 
1980) 

rii ~ n GJ min(l,C) , (2.15) 

with 

C = (ikev) (li) (WT ) (20) (se) (lO^kgrn" 3 ) 6XP (r~kff) • 

(2.16) 

The value of is very sensitive to the surface temperature for T s ~ T i: where Tj is the 
characteristic temperature 

where the cohesive energy of the condensed Fe 56 matter in a strong magnetic field is 
estimated from (2.6). A small change in T s around T s ~ Tj can have a large effect on 
rii, with a change by a factor of two causing ni/n GJ to vary from exponentially small to 
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approximately unity. For T s > Tj the existence of a magnetic metal at the neutron star 
surface does not affect the ejection of ions, and one has E\\ = on the stellar surface. 

Field emission of Fe 56 ions from a cold magnetic metal at the surface of pulsars with 
B s > 10 s T is unimportant for all known pulsars (Ginzburg and Usov 1972). Therefore, in 
the case T s < Tj there is no ejection of ions from the neutron star surface, and the electric 
field component E\\ near the polar caps is nonzero and is determined by the polar-gap 
structure (e.g., Ruderman and Sutherland 1975 and below). 



3. PARTICLE ACCELERATION AND CURVATURE RADIATION 

The acceleration of primary particles is due to the parallel electric field in the polar gaps. 
Curvature radiation is an important ingredient in the generation of the secondary pair 
plasma, and it can also limit the energy of the primary particles. 

3.1. Particle acceleration and curvature radiation 

If the energy losses of a particle are neglected, then a primary particle ejected at the stellar 
surface reaches a maximum energy 

e = Ymc 2 = eAip (3.1) 

at the top edge of the polar gap. Here V is the Lorentz factor and A(p is the potential 
across the polar gap. 

In a strong magnetic field near the pulsar surface, electrons lose the momentum com- 
ponent transverse to the magnetic field very rapidly (timescale <C R/c) and move away 
from the pulsar along the field lines. For such electrons in the ground-state Landau level, 
any energy loss is negligible up to the Lorentz factors of ~ 10. For 10 < V < 10 2 , the energy 
loss due to cyclotron resonant scattering of thermal X-rays from the neutron star surface 
increases sharply and may be very important for both motion of electrons and their radia- 
tion (Dermer 1990 and references therein). At very high energies, T 10 3 — 10 4 , the main 
energy loss for ultrarelativistic electrons in the pulsar magnetospheres is due to curvature 
radiation. In this case the rate of energy loss is (e.g., Ochelkov and Usov 1980) 

where R c is the radius of curvature of the magnetic field lines. Equation (3.1) is valid 
only if the rate of energy loss, is smaller than the rate of energy gain, e£j|, due to 
acceleration by the parallel electric field. The maximum energy of electrons cannot exceed 
the value implied by balancing these energy gains and losses, that is, the maximum is 
determined by |e| ~ e£j|C. The Lorentz factor of particles when this quasi-stationary state 
is achieved is 

r,r st =p^) 1/4 . (3.3) 
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Electrons in a strong magnetic field quickly radiate away their perpendicular momen- 
tum and are then in their lowest Landau level. Such electrons move along curved magnetic 
field lines and emit curvature radiation. The mean energy, f 7 , of curvature photons gen- 
erated by ultrarelativistic electrons, r>l, is 

For Aip = Aip ns , T = eA</? RS /mc 2 ~ 3 x 10 6 (Ruderman and Sutherland 1975) and 
R c ~ 10R ~ 10 5 m, from (3.4) we have e 1 ~ 10 2 MeV. Thus the curvature radiation of 
electrons accelerated in the polar gaps falls in the 7-ray range. Note also that the mean 
energy satisfies e 1 ^> 2mc 2 ~ 1 MeV. 

3.2. Curvature 7 rays and pair creation threshold 

A photon propagating in a strong magnetic field may decay into two photons or into pairs 
(Toll 1952; Klepikov 1954; Erber 1966; Adler 1971 and section 4 below). Pair creation 
occurs only when the energy of the photons exceeds the threshold ~ 2mc 2 in the frame in 
which the photon is propagating perpendicular to the field lines (section 4.2 below). 

Although a particle in its lowest Landau level has no perpendicular motion in a uniform 
magnetic field, it does have a small perpendicular motion if the field lines are curved. The 
perpendicular motion is at the curvature drift velocity, v c , which is such that the Lorentz 
force, oc v c x B, causes the particle to follow the curved field line (e.g., Zheleznyakov and 
Shaposhnikov 1979). The curvature drift implies a perpendicular momentum 

p cd = L_ , (3.5) 

mu) B R c 

where uj b = eB/mc is the gyration frequency. The angle, if), between the momentum 
vector p and the tangent to the magnetic field line, defined by p c d = psinifj, p\\ = p cos tp, 
V = IpI = r™ 7 is given by 

c;1.7x 10- (*) "f^m", (3.6) 
r u B R c V 108T / \l0 8 J\10 5 mJ ' v 1 

where (3.5) is used for ip <C 1. 

A highly relativistic particle emits nearly all its radiation in a cone with half angle 
~ r _1 <^ 1 about the momentum vector p. Hence, the curvature photons are emitted in a 
cone nearly parallel to the field lines, but centered on an angle ip away from the field lines. 
The question arises as to whether the condition for decay into pairs can be satisfied at 
the point of emission. The curvature photons have a perpendicular energy < £ 7 (^ + r _1 ). 
Using (3.4) and (3.6), this energy is below the threshold (2mc 2 ) for pair creation for 

3 cr 2 ( cr 2 \ 2mc 2 . n . 

+ 1 < ~z • (3.7) 



2 uj n R c \ uj b R c J hu) 
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This condition corresponds to the Lorentz factor of ultrarelativistic particles being smaller 
than 



r* = 



1 a ?\ 1/2 

1 4mc 2 \ 1 1 
4 3 huj B 



1/2 



(3.8) 



In the limits of a weak and strong magnetic field one gets 



AR^muj^ 

4R c mc 
3% 



1/4 



1/2 



for B < B CI , 
for B > B cr , 



(3.9) 



which gives numerical values 



R r 



10 5 



m 



1/2 



2.8 x 10 8 
5.8 x 10 s 



B 



10 8 T 



1/4 



for B < B cr , 
for S > S cr . 



(3.10) 



Inside the polar gap, the value of T is restricted by V < r st , where T st depends on E\\, 
cf. (3.3). The upper limit on E\\ when charged particles flow freely from the neutron star 
surface is (Arons 1981) 

and when there is no particle flow from the surface the upper limit is (Ruderman and 
Sutherland 1975) 



E 



RS 



VtR 



,) 
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(3.12) 



For all known pulsars one has E^ s ~ (10 2 - 10 4 ) Ep implying E^ s 3> ^j^. 
Substituting E^ s into (3.3) and using (3.9) we have 

for B s < S cr , 



r* Vo.is 




for B s > B 



(3.13) 



From (3.13) it follows that V is smaller than T* for all known pulsars except the millisecond 
pulsars. However, the magnetic field at the surface of millisecond pulsars, B s ~ 10 5 T, is 
not strong enough to determine the structure of matter, and the magnetic metal phase 
does not form in the surface layers of the neutron stars (section 2). Therefore, particles 
flow freely from the surface of such a pulsar. In this case, E^ is an upper limit on E\\ in 
the pulsar magnetospheres. Substituting E^ into (3.3) for En we get 



r r st 9 / P 

— < — < 2 x 10" 2 

r* r* Vo.is 



-5/S 



(3.14) 
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Then the condition P > 10 3 s, valid for all pulsars, in (3.14) implies V < 0.31T* for the 
millisecond pulsars. 

It follows that for all known pulsars the curvature photons generated near the neutron 
star surface are produced in a state below the pair creation threshold. In order to reach 
the threshold for decay into pairs (or the formation of positronium) the curvature photons 
must travel some distance so that the angle between their wave vector and the magnetic 
field has increased sufficiently. 

4. PROPAGATION OF 7 RAYS AND PAIR CREATION 

The propagation of 7 rays in a strong magnetic field can lead both to splitting of one photon 
into two or more photons and to decay of a photon into a free or bound electron-positron 
pair. These processes can depend on which wave mode the photon is in. 

4-1. The modes of propagation for 7 rays 

The conventional expression for the refractive index of a plasma, with the vacuum 
polarization by the magnetic field taken into account, differs from unity by the order of 
[0.1a(B / B CT ) 2 + (ujp/u) 2 ] sin 2 i?, with u p = (e 2 n p /eom) 1 ^ 2 , where n p is the plasma density, 
uj = e 1 /% is the photon frequency and d is the angle between the photon wave vector K and 
the magnetic field B (Erber 1966; Adler 1971; Shabad 1975). For u > 3a~ 1 / 2 (B CT / B)u pi 
the vacuum polarization gives the main contribution to the difference between the refractive 
index and unity. Near the pulsar surface, where B > 10 8 T and n p < 10 26 m -3 (Sturrock 
1971; Ruderman and Sutherland 1975; Arons 1979, 1981, 1983), this condition becomes 
u ^> 10 18 s _1 , which is well satisfied for 7-quanta. Hence, to understand the process of 
7-quantum propagation in the vicinity of a pulsar with B s > 10 8 T, it suffices to consider 
propagation in the vacuum polarized by a strong magnetic field. 

The principal modes of propagation for a photon in the magnetized vacuum are linearly 
polarized with electric vectors either perpendicular (_L mode) or parallel (|| mode) to the 
plane formed by the photon wave vector K and the vector B. The labeling convention 
adopted here is standard, although other labelings of the the modes are used (e.g., Adler 
1971; Usov and Shabad 1983). Both modes are generated in the process of the curvature 
radiation (Jackson 1975). 

Photons in either mode can decay into a free electron-positron pair provided that 
a relevant threshold is exceeded (Toll 1952). For the || mode both the electron and the 
positron can be in their ground states, and the relevant threshold is 

e 1 sin 1} = 2mc 2 . (4.1) 

For the _L mode the lowest allowed state requires that either the electron or the positron 
being in the first excited state, with the other in its ground state. The relevant threshold 
is then 

£ 7 sintf = mc 2 {l + [1 + (2B / B cr )] 1/2 } . (4.2) 
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4-2. The curvature photon splitting in a strong magnetic field 



While the photon is below the pair creation threshold, its main (inelastic) interaction with 
the magnetic field is a splitting into two photons 7+-B — > 7' '+7" +B (e.g., Adler et al. 1970; 
Bialynicka-Birula and Bialynicki-Birula 1970; Adler 1971; Stoneham 1979; Melrose 1983; 
Usov and Shabad 1983; Baring 1991 and references therein). For B < B cr , when radiative 
corrections are negligible, Bialynicka-Birula and Bialynicki-Birula (1970) showed that the 
splitting of a photon into more than two photons is unimportant. Adler (1971) showed 
that splitting involving an odd number of || -polarized photons (||— >|| + ||, ||— >-L + -L, 
_L— >|| + _L) is forbidden (by CP invariance), and that below the pair creation threshold, 
of the remaining transitions the only one that is kinematically allowed is _L— >|| + ||. The 
probability for _L-polarized photon splitting in the weak- field limit, B < £? cr , is greatest 
when the energy of the initial photons is divided equally between the two final photons. 

The optical depth for the splitting of _L-polarized photons with e 1 ^> mc 2 as they 
propagate from 1} ~ up to the pair creation threshold (4.2) is (Usov and Shabad 1983) 



r d ~ 1.7 x 10 £ 
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B rr 
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(4.3) 



It follows from (4.3) that if the strength of the magnetic field at the pulsar poles is high 
enough, B s > 10 9 T, most of the _L-polarized photons with e 1 < 10 2 MeV produced by 
curvature mechanism near the pulsar surface, are split and transformed into the || -polarized 
photons before the pair creation threshold is reached. However, the ||-polarized photons 
cannot be split below the pair creation threshold. 



4-3. Bound pair production in a strong magnetic field 

If the photon energy is above the pair creation threshold, the main process by which 
a photon interacts with the magnetic field is single-photon absorption, accompanied by 
pair creation: <y + B-+e + +e~+B (Klepikov 1954; Erber 1966; Tsai and Erber 1974; 
Melrose and Parle 1983; Shabad and Usov 1984). In the application to pulsars it is usually 
assumed that the curvature 7-quanta produced below the pair creation threshold propagate 
as polarized photons through the pulsar magnetosphere until they are absorbed by creating 
free pairs. However, before a curvature photon reaches the threshold for free pair creation 
it must cross the threshold for bound pair creation (Usov and Shabad 1985; Shabad and 
Usov 1985, 1986, cf. also Pavlov and Meszaros 1993 and references therein). Hence, when 
considering the creation of free pairs by curvature photons in the strong magnetic fields of 
pulsars, the process of conversion into bound pairs, that is, into positronium atoms, needs 
to be considered. 



4-3.1. Kinematic condition for positronium production 

To specify the positronium state created by a photon with a given momentum ftK, we 
choose a set of quantum numbers suggested by the special choice of gauge with scalar 
potential zero and vector potential with components A x = —By, A y = A z = for the 
constant and homogeneous magnetic field B = B z , B x = B y = 0. Since A does not 
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depend on x, z or t, the x- and ^-components of momentum and the energy are constants 
of the motion. Without loss of generality we choose the x-axis along ftKj_. Conservation 
of the ^-component of momentum implies that the parallel momentum of the center of 
mass of the positronium is equal to the parallel component of the photon momentum: 
TiK\\ = P z = p+ + p~ , where pf are the parallel momentum components of the electron 
and positron. 

For present purposes it suffices to consider the case B ^> a 2 B CT ~ 2.35 x 10 5 T, in 
which the orbital radius of the electron, r B = {%/ eB) 1 / 2 , is much less than the Bohr radius, 
a = AksqTi 2 /me 2 . In this case the dependence of the electron and positron wave functions 
on the perpendicular motion is of the same form as when they do not form positronium 
(Schiff and Snyder 1939). The perpendicular motions of the electron and positron are then 
each described by a continuous and a discrete quantum number. The continuous quantum 
numbers, p^, are interpreted as the y-coordinates of the gyrocenters y ± through 

= ("> 

Conservation of the x-component of momentum then implies hK± = P x = p+ + p~ ' . 
Thus, a photon with perpendicular momentum component HK± produces an electron and 
a positron whose gyrocenters are separated along the y-axis by HKj_/eB. The Coulomb 
interaction between the electron and the positron affects only the wave functions of the 
relative motion along the magnetic field, and reduces to the problem of the one-dimensional 
hydrogen atom (Loudon 1959). 

Thus, in a strong magnetic field the energy states of positronium in its rest frame may 
be labeled by five quantum numbers: n > and n' > for the Landau levels of the electron 
and positron, n c > for the principal quantum number of the one-dimensional hydrogen- 
like atom, the parity of the hydrogen-like atom state, and P 2 > for the separation of 
gyrocenters of the electron and positron. The energy of positronium does not depend on 
Px = Px~ ~ Pxi reflecting the fact that the energy cannot depend on the location of the 
center of mass of the positronium in a uniform magnetic field. 

The corrections to the dispersion relation for the two modes due to the polarization 
of the vacuum are not important in the present context. Hence we may approximate the 
dispersion relations, expressed as the energy of a photon, by the vacuum relation 

e 1 = hc(Kf + K 2 ± ) 1/2 . (4.5) 
The energy of positronium may be written 

e p =[P 2 c 2 +e 2 nn/ (n c ,P 2 h- 2 )} 

where 

e nn >(n c ,P 2 h~ 2 ) = mc 2 

is the rest energy of the positronium, with Ae nn >(n c , P 2 h~ 2 ) its binding energy. 

1 a 



(4.6) 



1 + 



2nB 



B r 



1/2 



+ 1 + 



2n'B 

Brr 



1/2 



- Ae nn ,(n c ,P 2 h- 2 ), (4.7) 



Using both momentum conservation, hK± = P x and hK\\ = P z , and energy conserva- 
tion, e 1 = s p , we obtain the kinematic condition for positronium production by a photon 
in the form of an equation for K±: 

K ± = ±-e nn ,(n c ,Ki). (4.8) 



In order to establish whether this equation admits any solutions, one should know the 
positronium dispersion law, that is, the energy dependence on P x . Unlike the case of free 
pair production by a photon in a magnetic field, to find the kinematic condition for bound 
pair production a dynamic problem has to be solved. This problem was considered by 
Shabad and Usov (1986) in detail, and the energy of a positronium atom was calculated 
for a state when all quantum numbers, n, n', n c , P z , P x and P x = p x — p~ , are arbitrary. 
Here the quantum number P£ = —eB(y + + y~) describes the location of the center of 
mass of the positronium in the y-direction, cf. (4.4). 



4-3.2. Energy of positronium and mixed photon-positronium states 

When the two photon dispersion curves (identical curves for each mode) and the energy 
relations for positronium for each of the discrete states are plotted in momentum space, the 
intersection points between the curves define parameters where the photon and positronium 
states interact. The interaction may be described in terms of the dispersion curves for the 
photons and the energy curves for the positronium states reconnecting to form a set of 
mixed states. Each of these mixed states is photon-like in one limit and positronium-like 
in another limit. A kinematic restriction implies that photons in the two different modes 
can interact only with a positronium state of the appropriate parity. Hence, there are 
actually two distinct sets of mixed states, each applying only to one photon mode and one 
parity state of the positronium. An adiabatic change, corresponding to evolution along 
the dispersion curve, then allows a photon in one mode to evolve into positronium in the 
appropriate parity state. 

Consider the simplest example, which is the ground state of the positronium, whose 
parity is positive and which interacts only with the || mode. For the ground state, n = 
n' = n c = 0, the binding energy is (Shabad and Usov 1986) 



Ae o(0, P 2 h~ 2 ) = a 2 mc 2 
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(4.9) 



From this equation, we can see that the binding energy is nonzero for any value of P x and 
tends asymptotically to zero as P 2 tends to infinity. Hence, if the magnetic field is strong 
enough, B ^> a 2 B cr , the positronium energy £oo varies in a narrow range from 

£ 00 (0, 0) = mc 2 {2 - « 2 [ln(a/r B )] 2 } , (4.10) 



to 2mc 2 when P 2 varies from zero to infinity. Therefore, (4.8) must have a solution. It 
follows that positronium production by a photon in the || mode becomes kinematically 
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allowed before the threshold of free pair production (4.1) is reached. A similar situation 
applies to a photon in the _L mode near the second threshold of pair creation (Shabad 
and Usov 1986). Figure 1 solves (4.8) graphically for both ||-polarized and _L-polarized 
photons. 

When the photon and positronium dispersion curves are treated as independent, they 
are given by the following equations (Figure 1): 




(4.11) 



for the photon, and 
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for the positronium with n = n' = n c = 0. In (4.12) the parallel and perpendicular wave 
vectors, K\\ = P z h~ x and K± = P x h~ x ^ are used instead of P z and P x to describe the 
positronium state. Mutual transformations of the photon and positronium are kinemati- 
cally allowed at the point satisfying K± = (K±)o, 

( ^ )0 .-{ 4 _ a2 ( ln _|_) 2 } 1/2 , (4 , 3 ) 



where, in the absence of any interference, the photon and positronium dispersion curves 
would intersect. 

The dispersion curves of the photon and the positronium interfere strongly near K\ = 
(K±)o- The dispersion curves repulse each other and reconnect to form mixed states 
(Figure 2). This effect is quite general (Von Neumann and Wigner 1964) and is sometimes 
called the theorem on nonintersection of spectral terms. 

For n = n' = n c = 0, the split dispersion curves of the mixed states are (Shabad and 
Usov 1985, 1986) 
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The lower of the two dispersion curves (the minus sign in (4.14)) is photon-like for K\ <^ 
(K±)q : and positronium-like for K\ 3> (K±)q (Figure 2). The upper branch in Figure 2 is 
positronium-like for K\ <^ (K±)q. This upper branch crosses the dispersion curve of the 
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odd state of the positronium with n c = 1 without interfering with it (Shabad and Usov 
1986). As K\ grows further, the upper branch approaches the even positronium state with 
n c = 1. There is an infinite set of bound positronium states at higher quantum numbers, 
and hence an infinite number of branches of the mixed-state dispersion curves below the 
limit, (e/hc) 2 — K 2 = (2mc/h) 2 , that separates the bound states from the continuum of 
unbound electron-positron states. 



4-3.3. Channelling along the magnetic field lines 

A curvature photon is emitted almost parallel to the magnetic field below the pair creation 
threshold, K\ <C (2mc/h) 2 (section 3.3). Provided that the geometrical optics or adiabatic 
approximation applies, the curvature of the magnetic field lines cause the photon to shift 
along the lower branch of the dispersion curves in Figure 2. As a result it passes through the 
mixed photon-positronium state and gradually turns into positronium when the condition 

K\ > (2mc 2 /h) 2 + AK]_ (4.16) 



is satisfied, where 
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is a measure of the gap between the lower and upper branches. 

In order for 7-quantum absorption to lead to the formation of bound pairs rather than 
free pairs, the adiabatic approximation must apply. This requires that the time (Shabad 
and Usov 1986) 
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needed to convert from the photon-dominated into positronium-dominated state should be 
much longer than the inverse frequency of the photon, u~ x = Uje 1 . For photon energies 
e 1 3> 2mc 2 which are of interest here, we have K 2 ~ (e 1 /%c) 2 . Taking this into account 

and using (4.15) - (4.17), the condition for the adiabatic approximation to apply, r c ^> cu _1 , 
can be written as a restriction on the magnetic field: 
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(4.19) 



For #0 ~ 10R ~ 10 5 m, (4.19) yields B > 0.03 B cr . 

Another restriction is imposed by a widening of the dispersion curves that occurs 
because of the finite lifetimes of all excited states. In the case of the || -polarized photons, 
for kinematic reasons the state described by the lower dispersion curve of equation (4.14) 
cannot decay into two photons (cf. section 4.2 and references therein). However, the 
positronium-dominated state which is described by the left part of the upper dispersion 
curve of equation (4.14), may decay into two photons. The probability of a two-photon 
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positronium annihilation in a strong magnetic field (B ~ 10 8 - 10 9 T) is (Wunner and 
Herold 1979) 



for the case K± = in the positronium rest frame. The value of W^^ decreases as K± 
increases, so that (4.20) gives the maximum value of W 2l - 

The widening of the dispersion curves can be neglected if the energy gap between 
the ± branches is much larger than this widening. From (4.14), in the rest frame of the 
photon-positronium state, K\\ = 0, the energy gap is (Shabad and Usov 1986) 



The widening may be neglected in evaluating UW^- The field which gives AsgjUW^ = 1 
may be found from (4.17), (4.20) and (4.21) to be about 0.08 B CT . The ratio Ae g /hW 2l 
is sensitive to the strength of the magnetic field: it is ~ 10 -4 for B = 0.05 B CT and ~ 10 
for B = 0.1B cr . Hence, for B > 0.15 cr the conditions for applicability of the dispersion 
curve (4.14) are fulfilled, and effectively all curvature 7-quanta in the || mode are captured 
near the first threshold of bound pair creation by gradually converting into positronium. 
Positronium atoms with n = n' = n c = and K\ > (K±)q, are stable against spontaneous 
decay into photons, which is forbidden by momentum and energy conservation (Herold et 
al. 1985; cf. also Adler 1971; Usov and Shabad 1983). Thus the positronium atoms formed 
by capture of the || -polarized photons in a strong magnetic field are stable, in the absence 
of such external factors as electric fields and ionizing radiation. 

The _L-polarized photons cross the first threshold (4.1) without any pair creation, and 
pair creation is allowed only when the threshold (4.2) is reached. Positronium created 
by _L-polarized photons at this threshold have n + n' > 1, so that either the electron or 
the positron is in an excited state with the spin quantum number, s = 1, opposite to 
that in the ground state. In this case the main reason for the widening of the dispersion 
curve is a spin-flip transition, s = l— > s = — 1. The probability of such a transition, W s , 
has been calculated for both the non-relativistic (Daugherty and Ventura 1978; Melrose 
and Zheleznyakov 1981) and the relativistic (Herold, Ruder and Wunner 1982) cases. The 
energy gap between the dispersion curves for a _L-polarized photon mixed with the positro- 
nium state of the second series, n = n c = 0, n' = 1 or n' = n c = 0, n = 1, differs little 
from (4.21). The condition for bound pairs to dominate the absorption of _L-polarized 
photons in a strong magnetic field is Ae g > hW s . Using the relativistic calculation of W s 
(cf. Figure 2 in Herold et al. 1982) and (4.21), we can write this condition as a restriction 
on the strength of the magnetic field: B > 0.15 B CT . 

Summarizing, the main results of this section follows. 

(1) For B < 0.04 B cr ~2x 10 8 T, the creation of free electron-positron pairs dominates. 

(2) For 0.04 B CT < B < 0.15 cr , the probabilities of creation of both positronium atoms 
and free pairs are more or less comparable and a numerical treatment is needed to 
determine the details. 

(3) For B > 0.1 - 0.15 B cr the curvature photons of both modes are captured near the 
pair creation thresholds. 




(4.20) 
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Note that when the condition B > 0.15 B cr is satisfied bound pair creation necessarily 
occurs. Bound pair creation occurs for lower values of £?, and dominates over free pair 
creation for B > 0.1 B CT . 



4-4- The role of strong electric fields 

So far in our discussion of bound pair creation we neglect the electric field in the pulsar 
magnetosphere. In general, both parallel, Ey, and perpendicular, Ej_, components are 
nonzero near the pulsar surface (Ruderman and Sutherland 1975; Arons 1979, 1981, 1983; 
Michel 1991). The field Ej_ causes a drift of particles across the magnetic field with velocity 

Vo = ^. (4.22) 

provided one has E± < B, which condition is satisfied in the cases of interest here. In the 
pulsar magnetosphere at a point with radius vector r, we may eliminate Ex by making a 
Lorentz transformation, with v = v D (r), to the drift frame. Quantities in the drift frame 
are denoted by a tilde, so that we have Ej_ = 0. 

The characteristic time, r D , for a test particle to be dragged into motion across B 
with the mean velocity v D is of the order of the time of one gyration: 

r D ~ u-'T ~ (10- 12 s) (-?—) 1 • (4.23) 
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For particles which are either ejected from the pulsar surface or created in the pulsar 
magnetosphere, the initial value of the Lorentz factor satisfies V < F st IS a few xlO 8 , 
cf. (3.3). As a consequence, the distance relativistic particles travel during the time r D 
is much smaller than the characteristic scale of the electric field variation, of order the 
neutron star radius. Hence, particles acquire the drift motion very rapidly. 

If the drift velocity v D is nonuniform its effect cannot be removed entirely by a Lorentz 
transformation. Suppose the drift velocity varies along the field lines, so that the electric 
field E± changes with distance, z. Then, in addition to the curvature drift (3.6), the po- 
larization drift contributes to the angle ip between the particle velocity v and the magnetic 
field B (Appendix A): 

^H^roy* ™ 

which applies to a particle in its lowest Landau level, and where Al is a characteristic scale 
length of the electric field variation. Emission of curvature photons is centered on the 
direction of the velocity of the particle, which is at angle ifj. For these curvature photons 
to be produced below the pair creation threshold requires 

e 7 sim/> < 2mc 2 . (4.25) 
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Assuming r ~ T, e 1 ~ e 7 , £? ~ £?, V < F st , E\\ < (-E| s ) max and E±/B < Oi?, and using 
(3.3), (3.4), (3.12) and (4.24), condition (4.25) gives 

9Tve ch f QR\ 5/2 R c < i 

v 2/5 

P>(1.6xlO- 3 s)^J . (4.27) 

For Al ~ R c near the neutron star surface, (4.27) is fulfilled for most pulsars with the only 
exceptions being a few millisecond pulsars. However, as mentioned above, the strength of 
the magnetic field at the surface of known millisecond pulsars is ~ 10 5 T, that is, much 
smaller than 0.1 B cr . For all known pulsars with strong magnetic fields, B s > 0.1 B cr , 
we have P ^> 10 -3 s, and hence the curvature photons generated near their surface are 
produced in a state below the pair creation threshold irrespective of existence of the field 
E± in the pulsar magnetosphere. Provided that the energy of the curvature photon is 
high enough, e 1 3> 2mc 2 , as the photon propagates it approaches the threshold for pair 
creation. With the characteristic scale of variation of the field E±_ is many orders of 
magnitude larger than the size of positronium atoms, at any instant we can introduce the 
frame with Ej_ = 0. In this frame our consideration of bound pair creation is applicable, 
implying that bound pairs must be created by the curvature photons for B > 0.1 B CI . 

The parallel electric field can lead to field ionization of the positronium. The Coulomb 
field and E\\ combine to form a potential barrier, and quantum mechanical tunneling 
through this barrier leads to a nonzero probability for the decay of the bound state into 
a free electron and a free positron (Herold et al. 1985; Shabad and Usov 1985; Bhatia et 
al. 1988). The probability of the field ionization in the rest frame of the positronium is 
(Shabad and Usov 1985) 
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The value of W E is very sensitive to the field E\\ . Provided that the strength of the field 
E\\ near the pulsar is less than 

the ionization probability of the positronium before it escapes from the pulsar environment 
is negligible. For E\\ > 2i?||° n , the mean free path of positronium atoms is much smaller 
than the neutron star radius, and positronium atoms are ionized almost immediately after 
their formation. Moreover, if the field E\\ is so strong that the characteristic time of the 
positronium ionization, (VFe) -1 , is of the order of or smaller than the characteristic time 
of the photon conversion, r c , into a bound pair, cf. (4.18) and (4.28), then free electron- 
positron pairs rather than bound pairs are created by the curvature photons. 

Numerically, E^ n varies slowly with the magnetic field and with the Lorentz factor of 

the positronium. For B s ~ (0.1 - l)B a and T ~ 10 - 10 6 , (4.29) gives Ef n ~ 10 12 Vm" 1 
to within a factor of 2 or so. 
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4-5. The photoionization of positronium atoms 

Positronium atoms can be destroyed through ionization by photons. The main source of 
such photons is thermal radiation from the neutron star surface. The surface temperature 
of neutron stars, with ages ~ 10 3 - 10 6 yr, is usually estimated to be ~ 10 5 - 10 6 K (e.g., 
Van Riper 1991; Page and App legate 1992 and references therein) These estimates are 
more or less in agreement with recent observations of X-rays from a few pulsars. The 
mean free path for relativistic positronium, T ^> 1, against photoionization is 

l ph ~cF(W ph )-\ (4.30) 

where W p h is the probability of the positronium photoionization in the rest frame. 

The value of W ph was calculated numerically by Bhatia et al. (1992) for B = 4.7 x 10 8 T 
and K± = (K±) ~ 2mc/%. The results of their calculations may be approximated by the 
analytic expression 

"^-< 6x1 ° 7s "'(if)" 2 (#k) 2 - < 431 > 

which is valid for 10 2 < V < 10 4 and 10 5 < T s < 10 6 K with the accuracy ~ 20 per cent. 
Equations (4.30) and (4.31) yield 

ip ^ (0 . 6xl03m) (-L) 3 (_^_) \ (4.32) 

Note that only photons with energy e 1 ~ (1 - 3)Aeoo in the rest frame contribute to the 
decay of the positronium (Bhatia et al. 1992). Hence, the dependence of W p \y on B and 
K± is through the B and K± dependence of the binding energy, which may be written in 
the form, cf. (4.9), 

a 2 

Ae OO {0, K 2 ± ,B) ~ — mc 2 

From (4.33) we infer that the B and K± dependence of W p h is only logarithmic, and hence 
(4.32) may be used to estimate the order of / p h for 0.1 B CT < B < B CT and K± > 2mc/h. 
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5. POLAR GAPS AND NONTHERMAL LUMINOSITIES OF PULSARS 

As discussed in section 4, adiabatic conversion of the curvature 7-quanta into mutually 
bound pairs in a strong magnetic field may result in these pairs not screening the electric 
field near the pulsar. The potential Aip across the polar gaps may then exceed A(p RS , 
leading to an increase in the theoretical estimate of the total energy flux carried by rel- 
ativistic particles from the polar gaps into the pulsar magnetosphere. In this section we 
explore how polar-gap models need to be modified to apply to pulsars with strong mag- 
netic fields to take account of the role played by bound pairs. Assuming free ejections from 



the stellar surface, we find that the modification gives at most a modest enhancement in 
the energy flux in primary particles, and this is inadequate to account for the observed 
nonthermal emission from some pulsars. Alternatively, assuming limited ejection from the 
stellar surface, we argue that there is a feedback mechanism that favors a self-consistent 
model. The energy flux in primary particles in this model appears capable of accounting 
for the observed nonthermal emission from most of the observed high-frequency pulsars. 



5.1. Free ejection of electrons from the polar cap 

First consider the case where electrons flow freely from the surface of the neutron star, 
which occurs when the temperature of the polar caps satisfies T s > T e , cf. (2.13). The 
model of Arons is then assumed to apply, with E\\ = at the stellar surface. In the 
outflowing plasma at distance r from the center of the neutron star, E\\ may be as high as 
(Arons and Scharlemann 1979) 

„ „ A . . / s/AR p at 0<s<AR p , 

£ »^» Smixi (i/rf/ 2 at s>AR p , ^ 

where % is the angle between the stellar magnetic moment and the spin axis, s = r — R 
is the distance from the surface of the star, AR P ~ (QR/c^'^R is the radius of the polar 
cap, AR P <C R, and is given by (3.11). Here and below we consider a pulsar with a 
dipole magnetic field. The potential which corresponds to the electric field E\\ is 

s 2 

r >/-5 l „ . ^ at s < AR V , 
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The total power carried away by both relativistic particles and radiation from the 
polar gap into the pulsar magnetosphere is 

L p ~ jVp rim eA<p , (5.3) 



where A^ pr i m is the flux of primary electrons from the polar cap, Acp = tp(R + H) — <p(R) is 
the potential across the polar gap and H is the thickness of the polar gap. Equation (5.3) 
is valid irrespective of whether the pairs created near the pulsar are free or bound. The 
version of pair creation determines only the thickness of the polar gap, H, and thereby 
affects the value of L p . 

Assuming that the pairs created by the curvature photons are free, the thickness of 
the polar gap, H = Hf, is (Arons 1981, 1983) 
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^^(O.OSs) (3^) (sinz) 2 / 7 (5.5) 

is the value of the pulsar period at which H is equal to AR P . 

The formation of bound rather than free pairs is important for B s > 0.15 cr . Near 
the polar cap of the pulsar the condition B > 0.1 B CT is satisfied up to a height 



Hu = R 



B« 



0.1 B r 



1/3 

- 1 



(5.6) 



above the pulsar surface. Thus, for Hi, > Hf, bound pairs rather than free pairs are 
created in the polar gaps. In this case, provided that the probability of ionization of the 
bound pairs near the polar caps is sufficiently small, the thickness of the polar gaps could 
satisfy H ~ H\,. It follows that for Hf, > Hf one needs to consider the ionization of bound 
pairs near the polar caps inside the layer < s < Ht when estimating the values of both 
H and L™^ x . 

For positronium atoms the probability of field ionization is negligible for E^ < E 1 ™ ~ 
10 12 Vm" 1 (section 4.4). Using (3.11) this condition may be written in the form 

2/5 

p>(a8xl(r2s) (ark;) ■ < 6 - 7 > 

which is satisfied for all known pulsars. 

The mean free path for photoionization increases rapidly with increasing Lorentz 
factor of the bound pair, Z ph oc T 3 , cf. (4.32). It follows that the photoionization of bound 
pairs with energy near the low-energy edge of their spectrum makes the main contribution 
to the density of free particles inside the polar gaps. This differs qualitatively from the 
case of B < 0.1 B CI for which the high-energy tail of curvature photons is responsible for 
both the creation of free pairs in the polar gaps and the screening of the electric field E\\ 
(Ruderman and Sutherland 1975; Arons 1981). 

Let us estimate the height, Hi, above the stellar surface at which photoionization of 
bound pairs is efficient enough to determine the structure of the polar gaps. The lower 
limit on the Lorentz factor of pairs created at the distance s from the pulsar surface is 
(Appendix B) 

4 / c \ V 2 R 



The energy of the curvature photons responsible for the creation of pairs with the Lorentz 
factor r m i n , hu = 2r m ; n mc 2 , is less than the mean energy of the curvature photons radiated 
by the primary particles in the polar gaps, hu <C (3hc/2R c )F 3 , where V is the Lorentz 
factor of primary particles accelerated in the polar gaps (section 3.1). In this case the 
spectral power of curvature radiation does not depend on T and is given by (e.g., Ochelkov 
and Usov 1980) 

2 / \ 2/3 
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At distance s from the pulsar surface the density of the curvature photons with energy 
hu = 2r m i n mc 2 is 

p(w) s 

T 



n^n GJ ^-. (5.10) 



Taking into account the motion of the pair, the probability of the photoionization of 
bound pairs in the frame of the pulsar is W p h/T. From (4.31) and (5.8) - (5.10) we have 
the following estimate of the density of free pairs, n±, created near the polar caps: 

n±~n 1 ^^-, (5.11) 
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where we assume R c = 10 5 m in (5.9). If at distance s the density of free pairs, n±, 
exceeds about n GJ (QR/c) (Arons 1981, 1983), the field E\\ is screened at distances beyond 
s. This condition together with (5.12) gives a lower limit on the height, H = i/™ m , for 
photoionization of bound pairs near the pulsar surface to be important: 
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In our estimate of the density of free pairs, n±, cf. (5.12), the energy 2r m ; n mc 2 is adopted 
for both the curvature photons and the created pairs. However, this energy is only the 
threshold for pair creation, and the actual mean energy may be a factor two or so higher. 
Alternative assumptions that give an estimate of the maximum value, H = i7™ ax , of Hi 
are that all curvature photons are generated by primary particles at s = 0, and all are 
emitted at threshold. If we assume that the mean energy of each pair is 4r m i n mc 2 , from 
(5.8) and (5.11) we have s ~ #™ ax ~ 3H™ in . We assume H, = 2H™ in . 

Our estimate of the thickness of the polar gaps is (for Hb > Hf) 

H= min{# 6 , max[H f ,Hi]}. (5.14) 



In summary, for free electron ejection from the polar cap, any increase in LJJj.fj* is 
rather small for all known pulsars, and it cannot be account for strong nonthermal X- 
rays and 7-rays observed from a few pulsars. The main reason for this is that E\\ in 
the outflowing plasma is rather weak, and the thickness, Hf, of the polar gaps is large 
even before taking into account its increase due to the creation of bound pairs. Hence, 
the potential Acp ~ E\\Hf is modest, and with at best a modest increase in H, it is not 
possible to increase A(p substantially. 



Oft 



5.2. Limited ejection of electrons 

Now consider the opposite case where electrons do not flow freely from the stellar surface. 
This requires that two conditions be satisfied simultaneously: thermionic emission of elec- 
trons must be ineffective, and ejection of electrons due to field emission must be ineffective. 
The former of these conditions requires T s < T e , cf. (2.13), and the latter requires E\\ < E e , 
cf. (2.14). Provided these conditions are satisfied, E\\ near the pulsar surface is of order 
the vacuum value E* s ~ (SIR/ 'c) 3 / 2 cB s , which is about 8y/3(c/tlR) ~ 10 2 - 10 4 times the 
maximum value of E\\ in the Arons model, cf. (3.11). 

Assuming that B is strong enough for the formation of bound pairs only to occur, a 
self-consistent solution for the polar gap may be constructed if field ionization of bound 
pairs is possible partially within the polar gap (see below). This requires E^ n < E^ s 

within the polar gap. The conditions E^ n < E^ s and E^ s < E e can be satisfied only 

for E [ ° n < E e . For B ~ (0.1 - l)B cr , comparison of (2.14) and (4.29) with (4.33) shows 

-E|l° n ~ 0.1 E e , implying that there can be a range, E]? n < E^ s < E e , where escape from 
the surface is ineffective and where field ionization of bound pairs is possible. 

When there is neither free ejection from the stellar surface nor ionization of bound 
pairs, vacuum conditions apply. Then the field E\\ near the polar caps should be (Ruderman 
and Sutherland 1975): 

*■<••'> -^K^hfe)*?*- (5 ' 16) 

where p is the distance to the magnetic axis. Granted that E\\ is unscreened, and so is 
given by (3.12), the condition E^ s = E 1 ™ defines a minimum period P = P± longer than 
which field ionization is unimportant: 

2^ / cBA 2/3 ^ 

In the following discussion we assume P < Pi, so that field ionization is important. 

There is an inconsistency if one assumes that E\\ is either much smaller or much larger 
than -E , || 0n - On the one hand, suppose that ejection of electrons from the stellar surface 
is severely limited, implying E\\ ~ E^ s ^> E^ n . For P < P\ bound pairs created in the 

region < s < ±AR p are ionized before they escape from this region. The resulting free 
positrons are accelerated toward the surface of the star, producing an intense flux that 
heats the polar cap. In the absence of ejection of electrons from the pulsar surface the 
temperature at the polar cap would be (Ruderman and Sutherland 1975) 



T ^^ 7K \oik) U) ■ (5 - 17) 

It follows from (2.13) and (5.17) that, for both B 3 = (0.1 - l)B cr and P < P u T RS is about 
an order of magnitude more than T e , as given by (2.13). This implies that thermionic 
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emission (due to the heating by positron bombardment) would produce an intense flux of 
electrons from the polar cap for temperatures well below T RS . This invalidates the original 
assumption that electron flow from the surface is severely limited. On the other hand, if 
one assumes E\\ <C E^ n for P < P\, then field ionization of bound pairs is negligible, and 
both the heating of the polar caps by reversed particles and the screening of the field E\\ 
are ineffective. As a result, one should have the vacuum field E\\ = E^ s , which invalidates 

the assumption E\\ <C E™ n . 

This inconsistency is avoided only if there is partial but not total screening of E\\ 
inside the polar gap. Specifically, this field needs to be restricted to E\\ ~ (1 - 2)E™ n . 
This is possible through a feedback mechanism. To see that there is a feedback, suppose 
E\\ decreases. This causes the rate of ionization of the bound pairs to decrease, so that 
the number of returning positrons heating the polar cap decreases, the temperature of 
the polar cap decreases and the rate of thermionic emission decreases. This leaves fewer 
(primary) electrons to screen E\\, which increases toward E^ s , providing the feedback. An 
analogous argument leads to the conclusion that a postulated increase in E\\ implies a 
sequence of processes causing E\\ to decrease. 

There is a self-consistent solution for the polar gap with a strong magnetic field, 
B(s = AR P ) > 0.1 S CT or Hf, > AR P (cf. Cheng and Ruderman 1980). The surface is 
heated to T s ~ T e , required for marginally effective thermionic emission. The number 
density of primary electrons, n pr i m , which is equal to n GJ in the limit E^ s ^> E^ n of 

complete ionization, is such that the difference n GJ — n pr i m is proportional to E^ s — E™ n . 

The potential difference is just E™ n times the height of the polar gap, which is identified 

as the Ruderman-Sutherland value \AR p . Thus this self- consistent model has 

( E \r\ 

T s ~ T e , n prim ~ n GJ ( l - -JL_ 1 , A<p ~ \E^AR, 
E\\ ~ G[s]G[±AR p - s]Ef n , 



(5.18) 



where Q[x] is the step function, which is equal to unity for x > and to zero for x < 0. 

To obtain (5.18) we assume that none of the parameters of the polar gap depend on 
time. Alternatively it may be possible to find a time-dependent gap solution, analogous 
to the sparking model of Ruderman and Sutherland (1975). We have not explored this 
possibility in detail. Our preliminarily considerations suggest that such a time-dependent 
model may imply a somewhat higher (than our time-independent model) total power 
carried away by both relativistic particles and radiation from the polar gap into the pulsar 
magnetosphere. 

Although photoionization of bound pairs occurs for P < Pi it has little effect on the 
model (5.18) provided the number density of ionized pairs, n±, at s ~ ^AR P is smaller 
than n GJ — n pr i m ~ n GJ (£ , ||° n /£ , |j iS ). From (5.12) and (5.18), the equality n± = n GJ — n pr i m 
defines a lower limit, P = P^, to the period: 
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For P < P 2 the number density of free pairs is sufficient to cause an additional screening, 
reducing E\\ to <C P|° n - Thus the model (5.18) applies only for P2 < P < P\. 

For P 2 < P < Pi, (5.3) and (5.18) imply that the total power carried away by both 
relativistic particles and radiation from the polar gap into the pulsar magnetosphere is 

L p ~ 7r(Ai2 p ) 2 n prim ceA V ^--^[l- j E rot , (5.20) 
where (e.g., Ostriker and Gunn 1969) 



is the rate of rotation energy loss of the neutron star. The maximum conceivable value 
of the energy flux in nonthermal particles is E TOt . This maximum is achieved (to within 
a factor of order unity) in the Sturrock model with L p ~ 7r(APp) 2 n pr i m ceA</? in (5.20) 
interpreted in terms of n prim = n GJ , cf. (2.1), and A<p = A</? max , cf. (1.1). However, as 
already remarked, the Sturrock model is not internally consistent. We now argue that the 
modified model (5.18) can produce L p close to this maximum. 



It is convenient to define the ratio rf = L p /E TOt of the spin-down power going into 



primary particles. In our modified model this fraction is given by 

r o ]j] lon ( E 1C 

^ = 41E_~£JL_ fi_J!_| at P 2 <P<Pi. (5.22) 



Equations (5.16) and (5.22) yield 
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(5.23) 



The value of rf^ has a maximum, rf 3 = 3/8, at P = 2 2 / 3 Pi ~ 0.6 Pi. 

For comparison with (5.22) or (5.23), suppose that both ejection of particles from the 
pulsar surface is limited and pairs are created free. Then the fraction of the spin-down 
power radiated by pulsars is (Ruderman and Sutherland 1975; Cheng and Ruderman 1980) 

/ R \ -8/ 7 / p \ 15 / 7 



From (5.22) and (5.24) we can see that for B s ~ (0.1 - l)P cr the ratio rj^/rj^ is ~ 20 at 
P = 0.6 Pi and ~ 10 2 at P = P 2 - Hence, at P 2 < P < Pi the nonthermal luminosity 
of pulsars with strong magnetic fields is considerably enhanced as a result of the creation 
of bound pairs instead of free pairs. The nonthermal luminosity of such a pulsar with 
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P ~ 0.6 Pi may be comparable with the spin-down power, as in the model of Sturrock 
(1971). 

The model (5.18) breaks down for P < P2 and for P > P\. Consider the lower 
limit, P = P 2 . According to (5.12) with (2.1), the density of free pairs resulting from 
field ionization increases with decreasing pulsar period, n± oc -P -5 , and for P < P 2 the 
number density of ionized pairs dominates in the screening of E\\, thereby reducing Aip in 
(5.18) and hence L p according to (5.19). The fraction of the spin-down power radiated by 
pulsars drops sharply from ~ 77^ for P > P 2 to ~ rjf for P <^ P 2 . For P > Pi and T s < T e , 
there is neither particle ejection from the pulsar surface nor a cascade of pair production 
in the polar gaps. In this case the nonthermal luminosity of pulsars is negligible. Hence, 
for pulsars with strong magnetic fields, the death line should be P = Pi (cf. Chen and 
Ruderman 1993). However, this predicted death line is not supported by observation, as 
discussed further in section 5.5 below. 

5.3. O • B < at the polar caps 

Ions tend to be extracted from the stellar surface in the polar caps of pulsars by any 
parallel electric field in the case O ■ B < 0. As for models with ft ■ B > 0, it is convenient 
to distinguish between the case where particles can flow freely from the surface due to 
thermionic emission, and the case where flow from the surface is limited. The discussion 
of these two cases closely parallels those in sections 5.1 and 5.2, respectively. 

Ions flow freely from the pulsar surface with density ~ n GJ if the temperature of 
the polar cap satisfies T 3 > T i: cf. (2.17). The potential, A(p, near the pulsar surface 
then cannot be much greater than Aip ns . The point is that for both Aip > Ap ns and 
H > Hi, cf. (5.13), the polar gap is unstable, in the sense that a cascade of pair production 
develops. The cascade is stopped only when either the potential becomes of order A<p RS 
or the thickness of the polar gap becomes ~ Hi when the photoionization of bound pairs 
is small. It follows from the discussion in section 5.1 that for H ~ H i: the potential across 
the polar gap is not much greater than Aip ns for any known pulsar. Any increase (due 
to the reduced screening resulting from the pairs being bound rather than free) in the 
energy flux in particles is at most modest. Hence, for free ejection of ions from the polar 
cap, the creation of bound pairs in a strong magnetic field cannot be responsible for the 
high-frequency radiation of pulsars. 

Provided that the temperature of the polar cap, T s , without heating by reversed elec- 
trons is smaller than Tj, thermionic emission of ions is negligible. The E\\ field distribution 
near the polar cap is then given by (5.15). As for the case described in section 5.2, for 
P < Pi, one has E\\ > E^ n and the potential Aip across the polar gap is substantially 
greater than the Ruderman-Sutherland limit, A<p RS ~ a few xlO 12 V. Such a gap is un- 
stable to pair production inside it. Indeed, if a pair is created inside the gap, the field 
E\\ accelerates the positron out of the gap and accelerates the electron toward the stellar 
surface. Since Aip > Ap ns , one (or both) of these particles is accelerated in the gap until 
it generates the curvature radiation which is absorbed inside the gap to create secondary 
pairs. The total number of such secondary pairs per primary particle is much greater than 
unity. With E\\ > E^ n these pairs are necessarily free. In turn, the secondary pairs are 
accelerated and generate curvature photons which are absorbed to create tertiary pairs 



and so on. This cascade ceases if the polar cap is heated by reversed electrons so that the 
surface temperature increases to ~ Tj. (One has Tj < T RS , cf. (5.17), and the heating can 
increase the surface temperature to T RS .) Then outflowing ions screen E\\, restricting it 
to E\\ ~ (1 - 2)E , ||° n . Our modified polar gap, summarized in (5.18), then applies, with 
T e replaced by Tj and with n pr im the sum of the density of outflowing ions and of the 
positrons created in the polar gap. Thus, (5.20) may be used to estimate the total power 
carried away by both relativistic particles and radiation from the polar gap into the pulsar 
magnetosphere. 

5.4- Interpretation of observational data 

Strong nonthermal radiation in X-rays and 7 rays has been observed from a few radio pul- 
sars (Gunji et al. 1994; Ulmer 1994). Some observational data and theoretical predictions 
on these pulsars are presented in Table 1. 

The strength of the magnetic field given in Table 1 is estimated from the slow down 
of the pulsar rotation (B s oc (PP) 1 / 2 ). This actually estimates the value of the magnetic 
field near the light cylinder (the inner boundary of the radiation zone). The strength of 
the magnetic field at the stellar surface depends on the assumption about the structure 
of the magnetic field, which is assumed dipolar. Thus, in Table 1, for each pulsar it is 
the dipole component of the magnetic field, B s , at the magnetic pole that is given. These 
values of B s are twice as strong as the estimates of Michel (1991), whose B- values apply 
to the magnetic equator, rather than to the magnetic poles which are more relevant for 
the polar gaps (cf. Shapiro and Teukolsky 1983). 

The values of B s given in Table 1 are lower limits on the actual surface fields. The 
estimated fields also depend on the angle, i, between the dipole axis and the rotation axis, 
and the quoted values for i = rr/2 correspond to the minimum value of the dipole moment. 
More importantly, the actual field may be much more complicated than a pure dipole. 
For example, suppose the magnetic field of a pulsar were purely quadrupolar rather than 
purely dipolar; then to explain the same rate of deceleration of the pulsar rotation, the 
surface magnetic field would need to be (c/QR) ~ 10 2 - 10 4 times stronger. For PSR1055- 
52 and Geminga, the dipolar components give fields smaller than 0.1 5 cr . If this were the 
actual field then pairs would be created free. In suggesting that bound pairs might cause 
an increase in the pulsar luminosity, we need to assume a surface field stronger than the 
dipolar value. The value B s ~ 6 x 10 8 T is indicated in parenthesis for these pulsars in 
Table 1. 

From Table 1 we can see that the following inequalities hold for all pulsars: <C 
and r]° hs < r(L. We conclude that the effect of adiabatic conversion of the curvature 7- 
quanta into mutually bound pairs in a strong magnetic field, B > 0.1 P cr , is able to 
increase the nonthermal luminosities of pulsars sufficiently to account for the observations. 
However, the condition P > P2, cf. (5.19), is not satisfied for two of the pulsars listed, 
specifically for the Crab and the Crab-like PSR 0540-69. Apart from these two, we con- 
clude that photoionization of bound pairs is unimportant in determining the polar-gap 
structure, so that (5.23) applies. Thus the nonthermal high-frequency radiation of these 
pulsars may be explained in terms of the modified polar-gap model which takes account 
of the creation of bound pairs. 
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The Crab-like pulsars (PSR 0531+21 and PSR 0540-69) have periods shorter than P^. 
These are also young pulsars, which may not have had time to cool sufficiently (to T s < T e 
or T 3 <Ti) for our model to apply. The high-frequency radiation from these pulsars cannot 
be explained in terms of our modified model. However, for these two Crab-like pulsars, 
the outer-gap model of Cheng et al. (1986a) seems satisfactory (e.g., Ulmer et al. 1994). 
Moreover, the 7-ray emission from the Crab-like pulsars may also be explained in terms 
of the slot gap model of Arons (1983). However, the slot gaps is an effective source of the 
energy for nonthermal radiation only for dipole-like magnetic fields. Otherwise, the total 
power carried away by both relativistic particles and radiation from the slot gap into the 
pulsar magnetosphere is suppressed by a factor of ~ (R c / R)(QR/c) 1 / 2 (Arons 1983). It 
is reasonable that the radius of curvature of the magnetic field lines near the surface is of 
order the stellar radius, R c ~ R, in which case for the Crab pulsar the total power from 
the slot gap is even smaller than the the total power from the polar gap in the standard 
non-modified polar-gap models (Ruderman and Sutherland 1975; Arons 1979, 1981; Cheng 
and Ruderman 1980; Mestel et al. 1985). Therefore, the outer-gap model of Cheng et al. 
(1986a) which is not so sensitive to the structure of the magnetic field near the pulsar is 
preferable as a model of the high-frequency radiation from Crab-like pulsars. 

In Table 1, to estimate the fraction, r/^, of the spin-down power, -Z? ro t, radiated by 
pulsars in the modified polar-gap model it is assumed that only one polar gap is a source of 
both high-energy particles and radiation in the pulsar magnetospheres. This assumption 
is reasonable for PSR 1706-44, PSR 1509-58 and PSR 1055-52 which have only one 7-ray 
pulse in their 7-ray light curves (e.g., Ulmer 1994). The other pulsars have two pulses 
in their 7-ray light curves. This indicates that in their magnetospheres both polar gaps 
operate as particle accelerators. In this case the actual values of 77^ are twice the values 
given in Table 1. This change strengthens the conclusion that the modified polar-gap 
model is able to explain the observed high-frequency luminosities of pulsars. 

The distance to Geminga is particularly uncertain and the value adopted (150 pc) is 
also indicated in parenthesis in Table 1. If the distance is less than ~ 50 pc, we have 
r/° bs ^ r;^, and the standard polar-gap models could describe Geminga's 7 ray and X-ray 
emission adequately (Harding, Ozernoy and Usov 1993). The high-frequency luminosity 
of Geminga may be explained in the modified polar-gap model if the distance to Geminga 
is up to ~ 250 pc. 

Recently, Sturner and Dermer (1994) and Daugherty and Harding (1994) have shown 
that double pulses with large interpulse separations seen in the Crab, Vela, and Geminga 
can be produced by hollow cone emission from a single pole if the rotation and magnetic 
axes are nearly aligned. In these models the solid angle of 7-ray emission is very small, and 
the standard polar gap models have no problem supplying the power necessary to give the 
observed fluxes of 7-rays. However, in this case the chance of observing any given pulsar 
from the Earth is not more than ~ 10 -2 . As noted by Daugherty and Harding (1994), 
in spite of the poor statistics, it seems hard to reconcile such a low probability of pulsar 
detection with either the fraction of 7-ray pulsars observed among the known supernova 
remnants, or the fraction which have a double-pulse structure. 

5.5. The death line at P = Pi 

The prediction that P = P\ should be the death line (no free pairs for P > P\) for 
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pulsars with B > 0.1 B CT is not supported by the data. The distribution of radio pulsars 
essentially ignores our death line. The radio emission is attributed to free pairs and hence 
it appears that there is some source of free pairs even for pulsars with B > 0.1B cr and 
P > P\. There are several possibilities within the framework of our model. One suggestion 
concerns the decay of positronium from the excited state formed from _L-polarized photons 
(section 4.3). This decay should produce some free pairs, due for example to two-photon 
decay (cf. Melrose and Kirk 1986) to the densely-packed states (n c — > oo) around zero 
binding energy for the ground state (n = n' = 0). Even a small branching ratio for 
generation of free pairs may suffice to account for the apparent lack of suppression of the 
radio emission. Two other possibilities for the formation of some free pairs are suggested 
by a critical examination of the assumption that the curvature photons emitted by the 
primary particles cannot decay directly into free pairs (section 3.2). In our discussion 
in section 3.2 we assume that all curvature photons are confined to a forward cone with 
half angle ~ 1/r. However, a small fraction, ~ 1/r, of the power is emitted outside 
this cone. Some of these photons should satisfy the condition for decay into free pairs at 
the point of emission, and these photons should be a source of some free pairs. Another 
possibility follows from the fact that for any radiation process that produces photons with 
energy above the threshold for pair creation, the analogous process in which the photon is 
replaced by a pair is allowed. Thus direct curvature emission of pairs is another possible 
source of free pairs. Photon-photon production of pairs is yet a further possibility. The 
efficacy of these processes needs to be explored, but is subject to the major uncertainty that 
the number of free pairs required to account for the radio emission is poorly constrained. 

We conclude that P = P\ should be regarded as a death line for high-frequency 
emission, but not necessarily for radio emission, for which the luminosity is many orders 
smaller than the spin-down power. However, our arguments on this point are speculative, 
and a quantitative analysis is desirable. 



6. CONCLUSIONS AND DISCUSSION 

In this paper we consider pulsars with strong magnetic fields, B s > 0.15 cr , at their 
surfaces such that bound rather than free pairs should form due to decay of 7 rays in polar 
gaps. In discussing various relevant aspects of the underlying physical processes, we make 
the following points: 

1) . The surface structure of magnetic neutron stars determines whether particles flow 

freely from the surface, leading to a model of the type discussed by Arons (1981, 
1983), or whether particles are strongly bound to the surface, leading to a model of 
the type discussed by Ruderman and Sutherland (1975) and Cheng and Ruderman 
(1980). For pulsars with O • B > 0, the requirement on the surface temperature, T s , 
for electrons to remain strongly bound is T s < T e , with T e given by (2.13). For pulsars 
with ft • B < 0, the requirement for ions to remain strongly bound is T s < T^, with 
Ti given by (2.7). 

2) . In polar-gap models the primary particles emit curvature photons which are initially 

at too small an angle to B to decay into pairs, even when the curvature drift and the 
polarization drift are included. 
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3) . For such photons emitted near the stellar surface, s = 0, the curvature of dipolar field 

lines allows decay into pairs with a minimum Lorentz factor r m ; n that decreases ocl/s 
with increasing height, cf. (5.8). 

4) . The photons can be either ||-polarized or _L-polarized. If the magnetic field is strong 

enough, B > 8 x 10 8 T, photon splitting tends to cause the _L-polarized photons to 
decay into || -polarized photons. 

5) . In a moderately strong magnetic field, B>4x 10 s T, photons evolve into a bound 

pair (positronium) rather than decaying into a free pair. The ||-polarized produce 
positronium in its ground state, and the _L-polarized produce positronium in an excited 
state that quickly decays to the ground state through gyromagnetic emission (actually 
a spin-flip transition). 

6) . The positronium atoms can be destroyed (a) by field ionization for E\\ > E l ^ n , with 

-E|l° n given by (4.29), or (b) by photoionization due to thermal radiation from the 
neutron star, with mean free path given by (4.32). 

In applying these properties to polar-gap models for pulsars we distinguish between two 
types of model mentioned in point 1). above: models in which particles flow freely from 
the surface, implying E\\ = at the surface (e.g., Arons 1981, 1983), and models in which 
particles are tightly bound to the surface, so that has the vacuum value just above the 
surface (e.g., Ruderman and Sutherland 1975). The parallel electric field in the former, E^ 
given by (3.11), is much smaller that that in the latter, E^ s given by (3.12), and there is a 
corresponding difference in the potential drop, Aip ~ E\\H, across the polar gap of height 
H. Our particular interest is in the energy flux in primary particles, which is assumed to 
determine the nonthermal, high-frequency luminosity (e.g., Harding and Daugherty 1993). 
The energy flux in primary particles is n^^cnAR^eAip, where n pr i m is the number density 
of primaries and AR p ~ (fi-R/c) 1 / 2 !? is the radius of the polar cap. For n pr i m ~ n GJ , cf. 
(2.1), and A(p ~ Aip max = ^E^ s AR p , cf. (1.1), this is of order the spin-down power of the 
pulsar. However, Arons- type models have a small E^ <C E^ s , Ruderman-Sutherland-type 
models have small H <C AR P , and both have Aip <C A</? max , so that the power in primary 
particles is severely limited. The formation of positronium, rather than free pairs tends to 
reduce the screening of E\\, and hence to increase H. 

We find the following specific results for models in which bound pairs are formed. 

7) . There is at most a modest increase in the power in primary particles in Arons- type 

models, and even after modification to include positronium formation, these models 
cannot account for the observed high-frequency emission. 

8) . In the Ruderman-Sutherland-type models, the formation of bound pairs can have a 

large effect on the model, greatly increasing the height of the polar gap. This model 
requires T s < T e or T s < Ti, cf. point 1). above, and then photoionization of the 
bound pairs is unimportant. 

9) . For P > Pi, cf. (5.16), corresponding to E^ s < E™ n , cf. point 6). above, there is no 

field ionization of the bound pairs. This suggests that P = Pi should be the death 
line for pulsars with B > 0.1 B CT . Observationally, radio pulsars seem to ignore this 
death line. In section 5.5 we suggest processes that should produce some free pairs 
even for P > Pi, and note that only a modest number of pairs is required to account 
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for the radio emission. 

10) . For P2 < P < P± we argue for a self-consistent, time-independent model in which field 

ionization generates just enough free pairs for the reverse flux of primary particles to 
heat the polar cap sufficiently to maintain thermionic emission at the rate required 
for partial screening to reduce E\\ from E* s to ~ -^f|° n - This model is summarized in 
(5.18). The modified model is only weakly dependent on whether electrons or ions 
are pulled from the stellar surface; for ions, replace T e in (5.18) by Tj. 

11) . For P < P2, cf. (5.19), the number density of pairs from the field ionization exceeds 

the number density of primaries, and with decreasing P the model rapidly approaches 
that without bound-pair formation. 

12) . In the modified model the ratio, 77^, of the power in primary particles to the spin-down 

power, cf. (5.23), can be close to unity. In this sense the model is similar to that of 
Sturrock (1971), but without the internal inconsistency in the Sturrock model. 

In summary, the most important modifications due to the formation of bound pairs rather 
than free pairs is for pulsars with strong surface magnetic fields, B 3 > 4 x 10 8 T, with cool 
surfaces, T s < 5 x 10 5 K, and periods in the range P2 < P < P\. The model proposed here 
could then account for the efficient high-frequency emission. 

The high-frequency pulsars are (Table 1) PSR 0531+21 (the Crab), PSR 0540-69, 
PSR 0833-45 (Vela), PSR 1706-44, PSR 1509-58, PSR 1055-52 and Geminga. Of these, 
the middle three satisfy all the requirements for our modified model, but the first two 
(Crab-like) pulsars have periods P < P2 and the last two have surface magnetic fields 
B 3 < 4 x 10 8 T. For the Crab-like pulsars, the photoionization of bound pairs inside the 
polar gaps should be strong, and the modified polar-gap model is then not applicable. The 
outer-gap model of (Cheng, Ho and Ruderman 1986a) seems a viable alternative for these 
Crab-like pulsars (e.g., Ulmer et al. 1994). We argue that the magnetic field is likely to 
have nondipolar components and these are not included in the estimate of B s in Table 1. 
With plausible values for a nondipolar component, the model may also apply to PSR 1055- 
52 and Geminga. Recently, Usov (1994) has shown that the outer-gap model of Cheng, 
Ho and Ruderman (1986b) is inconsistent with the available data on Geminga, implying 
a polar-gap model. For Geminga the actual distance is important in estimating whether 
or not the modified polar-gap model is viable, as discussed in section 5.4. Alternatively, 
Harding, Ozernoy and Usov (1993) argued that, provided the distance to Geminga is not 
more than ~ 40 - 50 pc, it is possible to explain both the X-rays and 7 rays in terms of 
the polar gap model of Arons (1979, 1981). 

We discuss some other details of the application to high-frequency pulsars elsewhere 
(Usov and Melrose 1994). We conclude by commenting on the polarization of the high- 
frequency radiation. According to 4). above, for B s >0.2£? cr , most of the _L-polarized 
photons with e 1 < 10 2 MeV produced by curvature mechanism near the pulsar surface, 
are split and transformed into || -polarized photons before the pair creation threshold is 
reached (section 4.2). As a result, the 7-ray emission recorded from the pulsar at energies 
£ 7 < 10 2 MeV should be linearly polarized. Near the maximum of the light curve, the 7- 
ray polarization may be as high as 100 per cent. The degree of 7-ray polarization should 
decrease toward the edges of the 7-ray pulses. This is because the emission zone recedes 
from the neutron star surface as the phase of 7-ray emission shifts from the maximum 



of the light curve, implying that the B value in the line of sight to the region of 7-ray 
generation decreases. By observing the polarization of the 7-ray emission of pulsars it 
would be possible to estimate the strength of the magnetic field near the pulsar surface. 



APPENDIX A: PITCH-ANGLES IN THE DRIFT FRAME 



To estimate the angle ip between the particle velocity v and the magnetic field B for 
a relativistic electrons in its lowest Landau level, let us first consider the nonrelativistic 
motion of electrons in crossed B and Ej_(t) fields. The polarization drift (e.g., Schmidt 
1966; Landau and Lifshitz 1971) implies a drift (e/muj 2 ,) dE±(t)/dt. For a relativistic 
electron in its lowest Landau orbital, the polarization drift may be obtained by making a 
Lorentz transformation from the rest frame. The momentum in the drift frame is 



where (dEj_/dt) c is to be evaluated in the frame in which the electron is at rest on the 
average. With (E±) c = TE± and (dz) c = T~ x dz (e.g., Landau and Lifshitz 1971), we have 




(Al) 




(A2) 



From (Al) and (A2), we obtain 




mT dE± 



(A3) 



- eB 2 dz 



where p ~ mcT is the particle momentum in the drift frame. 
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APPENDIX B: MINIMUM ENERGY OF CREATED PAIRS 

Let us consider a neutron star with a centered magnetic dipole moment, p,. In cylindrical 
coordinates z, p, cp, the components of B near the polar cap are 

B B B =0 mi) 

z ~ An z3' An z* ' v ' 1 J 

where /U is the magnetic moment, z-axis is directed along the magnetic axis, z is measured 
from the neutron star center, p is the distance to the magnetic axis, and <p is the azimuthal 
coordinate. The magnetic field lines are determined by integrating 

dz dp 

bT z = b- p > (B2) 

which with (Bl) and (B2) gives 

p = p (z/z f / \ (B3) 

where po/z^ 2 is a constant. (The exact result is p = (po/zQ^ 2 )(p 2 + z 2 ) 3 / 2 .) 

For a curvature photon which is emitted at the point with z = zq and p = po along 
the magnetic field, K || B, the angle between the wave vector K and the magnetic field B 
in the process of the photon propagation is 

Taking into account that po < (Oi^/c) 1 / 2 /? for open magnetic field lines near the polar 
caps of pulsars, z > R and z — R R, (B4) gives 

3 {QR\ 1/2 z-R 



*{z)<-[ — ) —. (B5) 



At a distance s from pulsar surface, with s ~ z — R, we have 

3 fVLR\ 1/2 s 



If the energy of curvature photons, e 7 , is smaller than 

e min = 2mc 2 /i?, (B7) 

the process of pair creation by these photons is kinematically forbidden. From (B6) and 
(B7) it follows that the Lorentz factor of particles which can be created at the distance s 
exceeds 

A / c \V2 R 
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Table 1: The properties of high-frequency pulsars with data taken from Ulmer (1994). The 
distances to PSRs 1706-44, from Taylor, Manchester and Lyne (1993), and to 1055-52, 
from Koribalski et al. (1995) are indicated in parenthesis. 



Figure Captions 



Figure 1: The dispersion curves for photon, (4.5), and positronium, (4.6), are shown in the 
absence of interaction (solid curves). The interaction is taken into account schematically in 
the dotted curves, for ||-polarized photons, and the dot-dashed curves, are for _L-polarized 
photons, near the relevant intersection points (circles). The positronium states are defined 
by (4.7) with P x = hKj_; the labels on the curves denote the quantum number n c and the 
parity under reflection along the z-axis (superscripts). [After Shabad and Usov 1986.] 

Figure 2: As in Figure 1, with the solid lines denoting the mixed (||-polarization/positive- 
parity) photon-positronium states. The dots denote higher order (in n c ) states that merge 
into a continuum near e = mc 2 . [After Shabad and Usov 1986.] 
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